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Abstract 



We consider the wave propagating in the Einstein & de Sitter spacetime. The covariant d'Alembert's 
operator in the Einstein & de Sitter spacetime belongs to the family of the non-Fuchsian partial differential 
operators. We introduce the initial value problem for this equation and give the explicit representation 
, formulas for the solutions. We also show the LP — L q estimates for solutions. 
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1 Introduction 



The current note is concerned with the wave propagating in the universe modeled by the cosmological models 
with expansion. We are motivated by the significant importance of the solutions of the partial differential 
equations arising in the cosmological problems for our understanding of the universe. While there exists 
extensive literature on the hyperbolic equations, the question of initial value problems for the wave equation 
in the curved spaces with singularities, and, in particular, in the Einstein & de Sitter spacetime, which are 
well posed and preserve many features of the classical waves, remains unresolved. 

The homogeneous and isotropic cosmological models possess highest symmetry, which makes them more 
amenable to rigorous study. Among them, FLRW (Friedmann-Lemaitre-Robertson- Walker) models are 
mentioned. The simplest class of cosmological models can be obtained if we assume that the metric of 
the slices of constant time is flat and that the spacetime metric can be written in the form 



ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ) , 



with an appropriate scale factor a(t). (See, e.g. ,[22].) The assumption that the universe is expanding leads 
to the positivity of the time derivative j^a{t). The time dependence of the function a(t) is determined by 
the Einstein field equations for gravity 

Rptv — -^g^R — SttGT^. 

The metric of the Einstein & de Sitter universe (EdeS universe) is a particular member of the Friedmann- 
Robertson- Walker metrics 



ds 2 = -dt 2 + a 2 {t) 



dr 
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where K — —1,0, or +1, for a hyperbolic, flat or spherical spatial geometry, respectively. The Einstein 
equations are simply 

„ , ,d a 4tt , „ . / a \ 2 8tt K _ . 

p=-3(p + p)~, - = - — (p + 3p), - = p - (i.i) 

a a 6 \a I A a z 



where p is the proper energy density and p is pressure. For pressureless matter distributions (p — 0) and 
vanishing spatial curvature (K = 0) in the EdeS universe the solution to the equations (1.1) is 

a(t) = a i 2/3 , 

where ao is an integration constant [19]. This model describes an open geometry (the K = and spatial 
sections are diffeomorphic to K 3 ) in the presence of a constant non-zero energy density distribution. Even 
though the EdeS spacetime is conformally flat, its causal structure is quite different from asymptotically flat 
geometries. In particular, and unlike Minkowski or Schwarzschild, the past particle horizons exist. The EdeS 
spacetime is a good approximation to the large scale structure of the universe during a matter dominated 
phase, when the averaged (over space and time) energy density evolves adiabatically and pressures are 
vanishingly small, as, e.g., immediately after inflation [35]. This justifies why such a metric is adopted to 
model the collapse of overdensity perturbations in the early matter dominated phase that followed inflation. 



The Einstein & de Sitter model of the universe is the simplest non-empty expanding model with the 
line-element 

ds 2 = -dt 2 + a 2 t 4 ^ (dx 2 + dy 2 + dz 2 ) 

in comoving coordinates [11]. It was first proposed jointly by Einstein & de Sitter (the EdeS model) [10]. 
The observations of the microwave radiation fit in with this model [9] . The result of this case also correctly 
describes the early epoch, even in a universe with curvature different from zero [5, Sec. 8.2]. Recently it was 
used in [29] to study cosmological black holes. The key observation for that approach is that the line-element 
can also be written in the conformally flat form 

ds 2 = T 4 [-dr 2 + dr 2 + r 2 (d9 2 + sin 2 6d<f> 2 )], 

where the timelike coordinates are related by dr/dt = (3t)~ 2 / 3 . The last form is an asymptotic for the 
Schwarzschild metric whose line element may be written in the form 



ds 2 = T 



^ d e + ^dtdr + (l + dr 2 + r 2 (de 2 + sin 2 9d0 2 ) 



where t and r are timelike and spacelike coordinates related to the standard Schwarzschild coordinates t and 
f, by 

t = t + 2m In 



2m 



The fact that the resulting metric is asymptotically Einstein & de Sitter with the source reducing to a 
comoving pure dust at null infinity is used in [29] . In this sense the solution could be interpreted as a black 
hole in the asymptotic background of the Einstein & de Sitter Universe. 

The covariant d'Alambert's operator in the Einstein & de Sitter spacetime is 
Consequently, the covariant wave equation with the source term / written in the coordinates is 

(I) 2 *— jc (£)"*+?!*=/■ m 

The last equation belongs to the family of the non-Fuchsian partial differential equations. There is very 
advanced theory of such equations (see, e.g., [18] and references therein). In this note we investigate the 
initial value problem for this equation and give the representation formulas for the solutions with any 
dimension n G N of the spatial variable x € W 1 . 

The equation (1.2) is strictly hyperbolic in the domain with t > 0. On the hypersurface t = its 
coefficients have singularities that make the study of the initial value problem difficult. Then, the speed of 
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propagation is equal to t 3 for every t G R \ {0}. The equation (1.2) is not Lorentz invariant, which brings 
additional difficulties. 

The classical works on the Tricomi and Gcllerstcdt equations (see, e.g, [4], [6], [8], [31]) appeal to the 
singular Cauchy problem for the Euler-Poisson-Darboux equation, 

A« = % + ^ f , ceC, (1.3) 

and to the Asgeirsson mean value theorem when handling a high-dimensional case. Here A is the Laplace 
operator on the flat metric, A := J2i=i ~§^?- 

We use the approach suggested in [32] and reduce the problem for equation (1.2) to the Cauchy problem 
for the free wave equation in Minkowski spacetime: w tt — Av = 0. To us, this approach seems to be more 
immediate than the one that uses the Euler-Poisson-Darboux equation. More precisely, in the present note 
we utilize the solution v — v(x, t; b) to the Cauchy problem 

v tt -Av = 0, t>0, iel", 

v(x,0) = tp(x,b), v t {x,0) = 0, [ ' 

with the parameter b G B C R. We denote that solution by v v — v ip (x,t;b). In the case of function 
(p independent of parameter, we skip b and simply write v v = v lp (x,t). There are well-known explicit 
representation formulas for the solution of the last problem. We write those formulas to make the present 
note self-contained. If n = 1, and ip(x,t) = f(x,t) G C°°(R x R), B = R, then 

v f (x,t;b) = ^{f(x + t,b) + f(x-t,b)}, fel, xgR. (1.5) 
For / G C°°(R™ x R) and for odd n = 2m + 1, to G N, 

v f (x,t;b) = -(j-) 2 — } / J(x + ty,b)dS y , (el, x G R", (1.6) 

where, = 1 • 3 • 5 • . . . • (n — 2), while for x G R n with even n = 2m, m G N, 

Vf ( x ,t;b) = -(--) 2 ^ Hx + ty,b)dV y , t G R, x G R™, (1.7) 

where, Cq™' = 1 • 3 • 5 • . . . • (n — 1). (See, e.g. Theorems 4.1,4.2 [24].) In particular, if / is independent of t, 
then Vf(x, t; b) does not depend on b and we briefly write Vf(x, t). 

The straightforward application of the formulas obtained in [32] to the Cauchy problem for equation 
(1.2) decidedly does not work, but it reveals a surprising link to the Einstein & de Sitter spacetime. To 
demostrate that link we note that the "principal part" of equation (1.2) belongs to the family of the Tricomi- 
type equations (in the case of odd I it is Gellerstedt equation): 

u tt - t l A u = , 

where I G N. According to [32] the solution to the Cauchy problem 

u tt - t l Am = f(x,t), u(x,Q) = ifo(x), u t (x, 0) = (fi(x), 
with the smooth functions /, ipo, and ipi, can be represented as follows: 

U (x,t) = ^^/(l-^r'vf^w^ (1.8) 



r 2 (7) 

rt M)-<t>(b) 

+ 2c k db drE(r,t;0,b)v f (x,r;b), x G R", t>0, 

Jo Jo 
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with the kernel 



Here F(7,7;l;£) is the hypergeometric function (see, e.g., [1]), while k :— 1/2, 4>{t) :— j^f, 7 := 2FP21 

k 1 

and Cfe = (k + l) _ fcTT2~' E + T . The equation with I G N, x G M, and / = is studied in [30] by means of the 
partial Fourier transform and the confluent hypergeometric function. That approach gives parametrix of the 
Cauchy problem and as a consequence, a complete description of the propagation of the C°°-singularities. 

Suppose now that we are looking for the simplest possible kernel E(r,t;0,b) (1.9) of the last integral 
transform. In the hierarchy of the hypergeometric functions the simplest one, that is different from the 
constant, is a linear function. That simplest function F (a, b;l;Q has the parameters a = b = — 1 and 
coincides with 1 + C- The parameter I leading to such function F (— 1, — 1; 1; () = 1 + ( is exactly the 
exponent I = —4/3 of the wave equation (and of the metric tensor) in the Einstein & de Sitter spacetime. 

It is evident that the first term of the representation (1.8), as it is written, is meaningless if 7 = —1. This 
indicates the fact that the Cauchy problem is not well-posed anymore for the equation with I = —4/3. The 
next theorem also shows how the "lower order term" of the equation (1.2) affects the Cauchy problem. The 
main result of this paper is the following theorem. 

Theorem 1.1 Assume that <pi G C^ +3 - l (R n ), i = 0,l, f(x, t) G C^ +2 (R n x (0, 00)), and that with some 
e > one has 

\d%f{x, t)\ + \td t dPf{x, t)\ < C a t £ - 2 for all x G W l , and for all small t > 0, 
and for every a, (3, \a\ < [^] + 2, \(3\ < [^} + 1. Then the solution ip = ip(x, t) to the problem 
j ilht - t- 4/3 A V + 2t~ l i> t = f(x, t), t > 0, x G K", 

is given by 



1/3 



ip(x,t) = -t 2 db dsbvfix^f^s-^fy^ + b 2 / 3 - s 2 ) (1.11) 

2 Jo Jo 

+ t^v^ {x, 3t 1/3 ) - 3t- 2/3 (d t v V0 ) {x, 3t 1/3 ) + lf v Vl (x, 3t 1/3 s)(l - s 2 )ds . 

2 Jo 

The theorem shows that one cannot anticipate the well-posedness in the Cauchy problem for the wave 
equation in the Einstein & dc Sitter spacetime. In fact, it gives a structure of singularity of the solution at 
the point t = 0, which hints at the proper initial conditions which have to be prescribed for the solution. 
The initial conditions prescribed in the previous theorem are the Cauchy conditions modified to the so- 
called weighted initial conditions in order to adjust them to the equation. For the Euler-Poisson-Darboux 
equation (1.3) one can find such weighted initial conditions, for instance, in books [4] and [25] as well as 
in the references therein. The existence and uniqueness of the solutions for the initial value problem with 
the weighted initial conditions for the Euler-Poisson-Darboux equation and for the equation (1.3) with the 
time-dependent c are proved in [7] by application of the Fourier transform in rr-variable, as well as some 
transformations which reduce the equation to the confluent hypergeometric equation. 

Theorem 1.1 can be used to obtain some important properties for the solutions of the wave equation in 
Einstein & de Sitter spacetime, which are inherited from the solutions of the wave equation in Minkowski 
spacetime. In particular, as a consequence of the previous theorem, in Section 3, for the initial value problem 
(1.10) with n > 2, / = 0, and <po = 0, we obtain the following L p — L q estimate 

||(-A)-'V(-,t)||L.(R») < CtK 2 -"(|-|))||^|| iP(Rn) , t>0, 
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provided that s > 0, 1< p < 2, 1 + 1 = 1, ±(n + l)(± - ±) < 2s < n(± - i), and n(± - ±) - 1 < 2s. Similar 
estimates hold for the problem with general (po and /. Thus, in the present paper we prepare all necessary 
tools that will allow us to study in the forthcoming paper the solvability of semilinear wave equation in the 
Einstein & de Sitter spacetime. Having in mind the scale invariance of the equation and also the applications 
(see, e.g., [16], [28]), special attention will be given to the self-similar solutions. Results analogous to those 
presented in this note have already proven to be a good tool in the study of self-similar solutions [35] . 

This note is organized as follows. In Section 2 we prove the main theorem and give some of its extensions 
(Theorems 2.1-2.2) that allow stronger singularity in the source term. Section 3 is devoted to the application 
of the main theorem, namely, to the derivation of the L p — L q estimates. 

The EdeS model recently became a focus of interest for an increasing number of authors. (See, e.g., [2], 
[11], [12], [13], [14], [15], [22], [23], [29] and references therein.) We believe that the initial value problem 
and the explicit representation formulas obtained in the present paper fill the gap in the existing literature 
on the wave equation in the EdeS spacetime. 

2 Proof of the main theorem 

If we denote 

C := d 2 t -i-S A+2t- 1 d t , S:=d?-t-^A, 
then we can easily check for t ^ the following operator identity 

t~ 1 oSot = C. (2.1) 

The last equation suggests a change of unknown function tjj with u such that tjj = t~ 1 u. Then the problem 
for u is as follows: 

' u tt - t-^ 3 A u = g(x, t), t>0, xG E™, 

< lim u(x,t) = ifo{x), x G E™, (2.2) 
lim (u t (x,t) + 3i~ 1/3 A ip (x)j = <Pi(x), x G WL n , 

where g(x, t) = tf(x, t). Therefore it is enough to find a representation of the solution of the last problem. 
We discuss it in three separate cases of: (f) with ip = fi = 0; (ipo) with / = and ipi = 0; (</?i) with 
/ = and ip = 0. 

We will use the following property of the resolving operator of the problem (1.4): if P = P(D X ) is 
a pseudo-differential operator and r = r(t) is a smooth function of time, then the operator r(t)P{D x ) 
"commutes" with the resolving operator of the problem (1.4). More precisely, the following identity can 
easily be verified: 

v TP(Dx)f (x,t;b)=T(b)P(D x )v f {x,t;b) for all / G C°°(R n x (0, oo)) . (2.3) 
The operator S belongs to the family of the Tricomi-type operators 

T := d?-t l A, 

where I G E. The Cauchy problem for such operators with positive I, that is for the case of multiple 
characteristics, is well developed (see, e.g., [33] and references therein). The fundamental solutions of the 
operator and the representation formulas for the solutions of the Cauchy problem are given in [32] . The 
results of [32] allow us to write an ansatz for the solutions of the equation of (2.2). This ansatz has been 
used in [17] to extend the range of admissible values of I to negative numbers for the problem with data 
on the hyperplane t ^ 0. Here we use this ansatz to consider the weighted initial value problem (2.2) with 
data on the plane t — 0, where coefficients and source term are singular. As we already emphasized, it is 
interesting that the case of I = —4/3, that is the case of Einstein & de Sitter spacetime, is an exceptional case 
in the sense that it simplifies the Gauss' hypergeometric function F{^, 7; 1; z) appearing in the fundamental 
solutions constructed in [32], to the linear function F(—l, — 1; 1; z) = 1 + z. 
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The case of (f). Assume that / e C^ +2 (M. n x (0, oo)) and for every given multiindexes a, /?, \a\ < [f ] + 2, 
I/? | < [f ] + 1, the following inequality holds 

\d%f(x,t)\ + \tdtd£f(x,t)\ <Ct e ~ 2 for all x e E", te(0,T), 
for small positive T. We have to prove that the solution to the problem 

■ £V = /, t > 0, x e E™, 



{ 



is given by 



lim = 0, lim Vt(a;,t) = 0, 



V»(a;,i) = 



t 2 / d& / ds6w / (x,3t 1/3 s;tt)(l + 6 2 / 3 - s 2 ) . (2.4) 

Jo Jo 

Here the function Vf(x, r; t) is given by (1.5), (1.6), and (1.7), if n = 1, n is odd, and n is even, respectively. 

It is sufficient to check the properties of the function u = u(x,t) = tip(x,t), which solves the equation 
Su = g with g(x,t) = tf(x,t). Hence, we can restrict ourselves to the representation 



u(x, t) 



2 1 '"' 



o Jo 



1-6 



1/3 



ds v g (x,(t){t)s;tb)(l + b 2/3 - s 2 ) 



1 /•! /■? 

= Y8j dl J drv g {x,r;l){4> 2 {t) + <t> 2 {l)-r 2 ) (2-5) 

and take into account the identity v g (x,r- 7 b) = bvf(x,r;b). Here <fi(t) := St 1 / 3 . First we prove that the 
integral is convergent and that it represents a C 2 (R n x (0, oo))-function. We will skip the subindex g in the 
remainder of the proof. It is evident that v(x, r; b) <G C 2 (R™ x (0, oo)) and that 

\d%v{x,r;b)\ + \d r v{x,r;b)\ + \d 2 v(x, r; b)\ + \bd b v{x,r;b)\ < Cb- 1+£ , 

if \a\ < 2. It follows 

\v{x,<f>{t)s;tb)\ < Ct- 1+£ b- 1+£ . 
Then we use the last inequality and the first formula of (2.5) in the following inequalities: 

,1 r-l-b 1 ' 3 

\u(x,t)\ < Ct 1+e b- 1+e db ds (1 + 6 2 / 3 - s 2 ) 

Jo Jo 

< C e t 1+£ . 

The first formula of (2.5) leads to the estimate for the derivative 
d 



(2.6) 



dt 



u(x, t) 



< 



'At 



dsv(x,(f>{t)s;tb)(l + b 2/3 - s 2 ) 



+ 



+ 



r-1 rl-b 
/ db 

JO Jo 

r-1 r-l-b^'° 

db dss(a r w)(a;,< ? !)(i)s;^)(l + 6 2 / 3 - s 2 ) 

Jo Jo 

fdbf 
Jo Jo 



dsb{d b v)(x,<j){t)s;tb)(l + 6 2/3 - s 2 ) 



that implies 



9 < * 

g-U(X,t) 



,1 r-l-b 1 ' 3 

/ b- 1+£ db / ds (1 + 6 2 / 3 - s 2 ) 
Jo Jo 

rl r-l-b 1 ' 3 

/ b- 1+e db 

Jo Jo 



+C e i 1 / 3+e / b- 1+£ db / d ss (l + 6 2 / 3 - s 2 ) 



,i ri-b 1 ' 3 
/ b- 1+£ db / ds (1 + b 2 ' 3 - s 2 ) . 
Jo Jo 



+C £ tf 



(2.7) 
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Thus, the estimates (2.6) and (2.7) lead to the initial conditions 

lim u(x, t) = 0, lim uAx, t) = 0. (2.8) 

It remains to verify the equation. For the derivative -§^u(x, t) we use the second formula of (2.5) and obtain 

§- t u(x,t) = 0'(t)l^%( a; >(t)-^);O(^ 2 W + 2 (O-(0W-^(O) 2 )^ 

i ft Mt)-4>{1) 
+2<j/(t)<f>(t) — J dl J drv(x,r;l) 

i , -./ ft i , w ft f<t>(t)-4>{l) 

= ig(0 2 (*)) J v(x,tf>(t)-<f,(l) ] l)<i>(l)dl + —{f(t)) J^dlJ^ drv(x,r;l). 



For the second order derivative ^u{x,t), since (</> 2 (£))'w(2;, 0; £)<^>(t) = g(x,t) we derive from the last 
equation 

d 2 , N 

+ ~k { ^ {t)) " j v (*>0(*)-#O;O0(O* + ^(0 2 (t))V(*) / v r (x,ct>(t)-cf>(l);im)dl 
+ I8 i4,2{t)) " J M J drv{x 1 r-l) + -{4, 2 {t))'cj ) '{t) J v(x, cj>(t) - <f>(l); I) dl . (2.9) 

By means of the second formula of (2.5) and equation of (1.4) for the function An we derive 

r<t>(t)-4>(l) 



Au(x,t) = — dl rfra r 2 «(x,r;/)(^ 2 (t) + ^ 2 (0-r 2 ) 

18 Jo Jo 



/o Jo 

It follows 

1 /•* 

\dl 



Au{x,t) = l<P(t) f d r v{x,<f>{t)-<j>{l);l)(j){l)t 

I ft i ft f <t>(t)-<t>(i) 

— — / a r u(a;,0;Z) (</> 2 (t) + ^ 2 (0) * + 77 / di / drrd r v(x,r;l). 
18 Jo 9 Jo Jo 

Since d r v(x, 0; Z) = 0, one more integration by parts yields 

Au(x,t) = \<t>{t) I d r v(x,<j>(t)-<j>(l);l)<j>(l)dl (2.10) 
9 Jo 

1 ft I ft r 4>(t)-4>(l) 

+- V {x,4>{t) - - <p{i))di - g y o a y o dr^^o. 

Hence according to (2.9) and (2.10) the application of the operator df — t _4 / 3 A to the function u = u(x,t) 
gives 

u tt {x,t) -t~ 4/3 Au(x,t) 

= g(x,t) 

J o v( x , 4>{t) - m iw) di + ^(^ 2 (t))>'(t) J Vr (x, 4>{t) - m i w) di 

+ 2-(^(t)y ( p'(t) J* v(x, m - 0(0; di 



= g{x,t) 
+ 



^ 2 {t))" J* v(x, <j>(t) - <f>{l)- l)4>{l) dl + l(0 2 (i))>'(t) £ v(x, cj>(t) - 0(0; dl 
-t~ 4/3 l f Q <M*) ~ 0(0; O(0(*) - 0(0) <fi 

= g(x,t) + ±{4> 2 (t))'<l>'{t)J v(x,<P(t)-cP(l);l)dl-t- 4 ^cP(t) ^ v(x,cj>(t)-cj>(l);l)dl 
= 9(x,t). 

Thus, for this case the theorem is proven. 

One can allow more strong singularity of the source function / at t = 0. More precisely, one can reduce 
the case with such singularity to the one of Theorem 1.1 if the initial condition are modified. That is done 
in the next theorem. 

Theorem 2.1 Assume that f(x,t) G C^J+ 4 (R n x (0,oo)), t 2 f(x,t) G C(R n x [0,oo)), and that 

\d£f(x,t)\ + \td t dPf(x,t)\ < C a t- 2 for all t G (0,T), x G R™, 

and /or eijery a, /?, |a| < [§] + 4, |/3| < [§] + 3. Denote fo(x) := lim t 2 /(x, i) and suppose that with some 
e > /or i/ie functions f(x,t) and /o(x) £ C^l +4 (R n ) the following inequality is fulfilled: 

\d% (t/(a;,t) - t^foix)) | + |t9 t af (tf(x,t) - t^fax)) \ < C^ 1 for all t G (0,T), x G R™, 

/or every a, [3, \a\ < [f ] + 2, \(3\ < [f ] + 1. 
T/ien f/ie solution ip = ip(x, t) of the problem 



Vtt - i~ 4/3 A V + 2t- 1 ^t = /(&, t), i>0,i£ R", 

lim (tip(x, t)) = 0, Iim(*^ t (x,t)+V'(a;,i)-/o(ar)Ini)=0, xG 
t->o t-s-o 



(2.11) 



is given &?/ 



1 1 



t 18* 



t r <p(t)-Hi) 



o JO 



^(s,t) - -h{x)T{t) + — dl dr(0 2 (t)+0 A (l)-r 2 ) (2.12) 



x (z V/ (x, r; - r (x, r) + r 4 / 3 r(0 A w /o (x, r)) 

where r(t) := J Q * In sds. 

Proof. Consider the new unknown function w(x,t) := u — /o(x)r(t). Then 

$«;(*, i) - tf(x,t)-S(f (x)T(t)) 

= tf(x,t) - (t-'foix) - t- A l\(t) A /o(x)) 

= MM), 

where we have denoted 

Mx, i) := t (/(x, i) - t- 2 /o(x)) + t- A l\{t) A /o(x) . 
According to the condition of the theorem with some e > we have 

\d£h(x,t)\ + \tdtdgh(x,t)\ < C a t e - Y for all t G (0,T), x G M™, 
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a , Pi M < [f ] + 2, \[3\ < [^] + 1, that allows us to write representation (2.5) for the solution w — w(x,t): 

i ft f 4>(t)-4>(i) 

w(x,t) = — dl drv h ( Xl r-l){4> 2 (t) + 4> 2 (l)-r 2 ). (2.13) 

18 JO JO 



On the other hand, according to Theorem 1.1, the function w — w(x,t) satisfies initial conditions 
Consequently, 



limw(x,t) = 0, lim w t (x,t) = 0. (2.14) 
t^o t-X) 



lim u(x,t) = lim(w(x, t) + fo(x)r(t)) = 0, 
lim(u t (x, t) — fo(x) hit) = lim w t (x, t) = . 

For the function ip = tp(x, t) = t~ 1 u(x, t) this implies the initial conditions of (2.11). To prove representation 
formula (2.12), we note that 

v h (x,r;b) = v t{f{Xit) _ t -2 fo{x))+t -A/3 T(t)Mo{x) (x,r\b) 

= v tf(X)t) (x,r;b) - v t -i fa(x) (x, r; b) + v t -4/3 T{t)Afo{x) {x, r; b) 

= bvf(x, r; b) — b~ x vf a (x, r; b) + b~ 4 / 3 r(b) A Vf (x, r; b) 

= bvf(x, r; b) — b~ 1 Vf (x, r) + 6~ 4 / 3 t(&) A Vf (x, r) . 



Then we use (2.13) to write 

w(x,t) = —J dlj dr(cb 2 (t) + cf, 2 (l)-r 2 ) 



'0 Jo 

X 



(lv f (x, r; I) - l^Vf, (x, r) + r 4 / 3 r(0 A v fo (x, r)) . 
Thus, the representation 

u(x,t) = fo(x)r(t) + - j dl j dr(0 2 (t) + ( f, 2 (l)-r 2 ) 

x (lv f (x, r; I) - l^Vf^x, r) + r 4/3 r{l) A v fo {x, r)j 

implies (2.12). Theorem is proven. □ 

The last theorem does not exhaust the possible singularities of the source terms. The next theorem gives 
behavior of the solution as t — > if the source term is more singular. 

Theorem 2.2 Assume that f(x,t) e C^ +4 (R n x (0,oo)) and that with number a G (2,8/3) one has 
t a f{x,t) 6 C{W l x [0,oo)) and 

\dgf(x,t)\ + \tdtd£f(x,t)\ < C a r a for all t G (0,T), x G K", 

and for every a, [3, \a\ < [§] + 4 7 \(3\ < [^] + 1. Denote fo(x) := limt a /(x,t) and suppose that with some 

e > for the functions f = f(x,t) and fo = fo(x) G C^l +4 (R n ) the following inequality is fulfilled: 

|c£ (tf(x,t) - t x - a f {x)) | + \tdtd% (tf(x,t) - t^Mx)) | < Cat 6 " 1 for all t 6 (0,T), x G K", 

and for every a, /3, \a\ < [§] + 2, \(3\ < [§] + 1. Denote r{t) := (3 - a)- 1 {2 - a)-H 3 - a . 
Then the solution ip = ip(x,t) of the problem 

' tp tt - t- A ' 3 Aip + 2t- 1 xp t = f(x, t), t > 0, x G E", 
lim(^0- (3 _ a ; (2 _ a) ^/oW)=0, *€R», (2J5) 

lim ( tVt(x, t) + t) - ^— t 2 - a f (x) ) = , iel", 
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is given by 

1 i /-t r <t>(t)-<t>(i) 

1>{x,t) = -h(x)T{t) + j^J dl J dr(cj> 2 (t) + 4> 2 (l) - r 2 ) (2.16) 

x (lv f (x, r; I) - l 1 -^^ (x, r) + 1~ a I 3 t{1) A v fo {x, r)) . 

Proof. Consider the new unknown function w(x,t) := u — fo(x)r(t). Then 

Sw(x,t) = tf{x,t)-S{fo(x)T{t)) 

= tf(x, t) - ft-'Mx) - r 4 / 3 r(t) A /„(*)) 
= M^*), 

where we have denoted 

h(x, t) := i (/(a;, i) - t~ a f (x)) + t~ A '\{t) A f (x) . 
According to the condition of the theorem with some e > we have 

\d*h(x,t)\ + \tdght(x,t)\ < Ccf- 1 for all t € (0,T), ieP, 

a , /8j |«| < [§] +2, |/3| < [^} + 1, that allows us to write representation (2.13). On the other hand, according 
to Theorem 1.1, the function w — w(x,t) satisfies initial conditions (2.14). Consequently, 

\\va(u(x,t) — f (x)T(t)) = limw(x,t) =0, 

lim (u t (x,t) - f (x)T'(t)) = lim w t (x,i) = 0. 

t ■ 7 U T r U 

For the function = ij){x,t) = t~ 1 u(x,t) this implies the initial conditions 

]im(tip(x,t) - fo{x)r(t)) = 0, lim {pl> t (x, t) + ip(x, t) - f ( x y(t)) = , 

which coincide with ones of (2.15). To prove representation formula (2.16), we note that 

v h (x,r;b) = v t(f{xt) ^ t -a fo{x))+t -4/3 T(t)Afo(x) (x,r:b) 

= v tf{Xtt) {x,r;b) - v t i- afo{x) (x,r;b) + v t -4/3 T{t)Afo{x) {x,r;b) 
= bvf(x,r;b) - b 1 ~ a Vf a (x, r) + 6~ 4/3 t(6) A Vf (x,r) . 



Then 



1 



w(x,t) = ^ / dl I drv h ( Xl r-l){4> 2 {t)+4> 2 (l)-r 2 ) 



4>(t)-<j>(l) 

dr (0 2 (t) + 2 (Z)-r 2 ) 



x (lv f {x,r-l)~l 1 - a v ta {x,r)+r i l\{l) Av fo (x,rfj . 

Thus, the following representation 

i rt M)-4>{1) 
u(x,t) = fo(x)r(t) + — J dl J dr(0 2 (t) + 2 (O-r 2 ) 

x (lvf{x, r; /) - p-'vfrix, r) + r 4/3 r(l) A v fo {x, r)) 

for the function u = u(x,t) implies (2.16). Theorem is proven. □ 

The case of (fo). In this case / = and ip\ — 0. One can find in the literature different approaches for 
the construction of the solutions of the Fuchsian and non-Fuchsian partial differential equations. (See, e.g. 
[18], [20].) The next two lemmas give for / — behavior of the solutions of the equation of (1.10) near the 
point of singularity t = of the coefficients. 
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Lemma 2.3 For (p £ C [ * ]+3 (R n ) the function 

u(x,t) = v V0 (x,3t^ 3 ) -3t 1 / 3 (9 r t; V0 )(a;,3t 1 / 3 ) (2.17) 

solves the problem 

( Su = 0, x £ R", i > 0, 

1 limu(x,t) = <p (x), lim (u t (x, t) + St^ 1 ^ 3 A <po(x) \ = 0, iel" 

Here v v (x, it 1 / 3 ) is the value of the solution v(x, r) to the Cauchy problem for the wave equation, v rr — Av = 0, 
v(x, 0) = tp(x), v t (x,0) — 0, taken at the point (x,r) = (x^t 1 ^ 3 ). 

Proof. We verify it by straightforward calculations. It is evident that 

Au(x,t) = Aw V0 (x,3t 1/3 ) -3t 1/3 (9 r Av V0 (x,r)) r=3tl/3 . (2.18) 

Denote 

v (x,t) = v V0 (x,3t 1/3 ), vi{x,t) = -3t 1/3 {d r v Vo (x,r)) r=3tl/3 . 
Then, for the derivatives d t Vo(x,t) and d 2 vo(x,t) we have 

= f- 2 /3. 



d t v (x,t) = t /3 (d r v Vo (x,r)) r=3tl/3 , 
2 



d 2 v (x,t) = -\r 5 ' 3 (d r v m (x,r)) r=3tl/3 +r i ' 3 (d 2 r v V0 (x,r)) r=3tl/3 



At the mean time for the derivatives dtV\(x, t) and d 2 vi(x, t) we have 

d t vi{x,t) = -tr 2/3 {d r v Vo (x,r)) r=3tl/3 - 3i~ 1/3 (d 2 v Vo (x,r)) r=3tl/3 > 
2 

df Vl (x,t) = -t~ 5/3 (d r v V0 (x,r)) r=3tl/3 -3*" 1 (d 3 v Vo (x, r)) r=3tV , . 

Hence, for the first order derivative d t u(x,t) and for the second order derivative d 2 u(x,t) we have 
d t u(x, t) = -?,r 1/3 (d 2 v V0 (x, r)) r=3tl/ 3 , 

d 2 u(x,t) = f- 4 / 3 (a r V(a ; ,r)) r= 3 tl/ 3-3i- 1 (a 3 ^o(^.^)) r= 3 t i/3, (2.19) 
respectively. Consequently, using (2.18), (2.19), and the definition of v v we obtain 

d 2 u(x,t) - t-^ 3 Au( Xl t) = r 4 ' 3 (d 2 v VQ (x,r)) r=3tl/3 - M- 1 (d 3 v Vo (x,r)) r=3tl/3 

-t~ 4/3 {Av V0 (x,r)) r=3tl/3 +3t- 1 (d r Av V0 (x,r)) r=3tl/3 

= r i/3 (d 2 v V0 (x,r) - Av Vo (x,r)) r=3tl/3 

-St- 1 (d r (d 2 v V0 (x,r) - Av VQ {x,r))) r=3tl/3 

= 0. 

Thus, the function u = u(x,t) solves the equation u u (x, t) — i~ 4 / 3 Au(x,t) = 0. Lemma is proven. □ 
Corollary 2.4 The function ip — t- 1 u(x, t) solves the problem (1.10) with ipi = and with f = 0, that is 
j ij u {x, t) - t- 4 / 3 A $(x, t) + 2t- 1 ip t {x, t)=0, 

I ]imtip(x,t) = tp , lim (tip t (x,t) +ip(x,t) +3t~^ A <p (xj) = 0, x £ R™. 

\ t-s-0 t->0 V / 

In particular, the corollary shows that for the given dimension n £ N Huy gens' principle is valid for some 
particular waves propagating in the Einstein & dc Sitter model of the universe if and only if it is valid for 
the waves propagating in Minkowski spacetime (cf. with [26], [32], [36]). 

The case of (<£i). In this case / = and (po = 0. 
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Lemma 2.5 For <pi G C^ ]+2 (R n ) the function 



u(x, t) 



3 f 1 
^ JO 



(2.20) 



solves the problem 



(Su = 0, i>0, 

| lim i) = 0, lim Ut(x, t) = tpi(x), 



x e 



Here v v (x, (j>(t)s) is the value of the solution v(x, r) to the Cauchy problem for the wave equation, v rr — Av = 
0, v(x, 0) = (fi(x), v t (x,0) = 0, taken at the point (x,r) = (x,(j)(t)s), while <j>(t) = 3t 1//3 . 

Proof. We prove the lemma by straightforward calculations. We have 

3 f 1 1 

u{x,t) = t- J v Vl (x,4>(t)s)(l - s 2 )ds = — J v Vl (x,r)((t) 2 (t) - r 2 )dr. 



For the first order derivative we derive 

r<K*) 



dtu(x. 



I r<t>{t) i r4>(t) 

t) = d t — J v Vl (x,r)(<f> 2 {t)-r 2 )dr=-t- 1 ^ J v Vl {x,r)dr , 



while for the second order derivative using the last equation and integration by parts we obtain 



d 2 u{x,t) 



1 f n ' 1 

= -Q f ~ 4/3 I v Vl (x,r)dr + -t- 1 v ipi (x,<l>{t)) 

1 / I m f+ {t) \ 1 

= ~gt 4/3 I v Vl (x,r)r -J r(d r v ipi )(x,r)dr I + -t 1 v, 

1 , ( r m \ 1 

= -g*" 4/3 Ki(*>0(*M*) -J r(d r v Vl )(x,r)dr\ + -t~\ Vl (x, <j>(t)) . 



(x,4>(t)) 



Consequently, 



d?u(x,i) 



i r<Kt) 
= 9*" 4/3 y o r(d r v Vl )(x,r)dr. 



At the same time, we have 



1 r4>(t) 

Au(x,t) = — j o Av Vl (x,r)(4> 2 (t)-r 2 )dr. 



Then equations (2.21) and (2.22) imply 

u tt (x,t) -i~ 4/3 Au(x,t) 



- _y--4/3 
" 9 



{d r v Vl ) (x,r) 



i r4>(t) 

t^- J Av Vl (x,r)^ 2 (t)-r 2 )dr 



4>(t) 



{d 2 v Vl ) (x,r)dr 



1 -4/s U 2 (t) 
~ 9 2 



(<9 r v Vl ) (a;,0(t)) - y — (d 2 v Vl ) (x,r)d. 



(2.21) 



(2.22) 
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i i r<t>(t) 

= -^t-^4> 2 {t) (d r v Vl )(x,m — t- 4 / 3 ^ r 2 (%v Vl ) (x,r)d 



1 







= -t- 2 ^(d r v Vl )(x,cP(t))- z t ~ 2/3 J o A v Vl (x,r)dr . 
The definition of the function v ipi suggests that the function u — u(x,t) solves the equation: 

i ( r m 1 

u tt (x,t) -t- 4 / 3 Au{x,t) = -r" 1 ' 3 I (d r v Vl )(x,<l>(t))- J {d 2 r v Vl ){ Xl r)dr\ 

= 0. 

Finally we verify the second initial condition by means of the l'Hospital's rule: 

1 _ 

limut{x,t) = lim -t 1/3 / v Vl (x, r)dr = Yimv Vl {x,(j)(t)) = v Vl (x, 0) = (pi (x) . 

t — >-0 t — >-0 o J q t — 

Lemma is proven. □ 

Corollary 2.6 The function = t~ 1 u(x,t) solves the problem (1.10) with ipo = and without source term 
f = 0, that is 

j ^ tt - i-4/3 A V + 2*-^ = 0, O0,i6R", 

] lim tip(x, t) = 0, lim (tip t (x,t) + ip(x,t)) = (fi{x), x 6 R n . 

The last corollary completes the proof of Theorem 1.1. □ 
In particular, Corollary 2.4 and Corollary 2.6 show that, because of the integration in the formula (2.20), 
for all n G N Huygens' principle is not valid for waves propagating in the Einstein & de Sitter model of the 
universe, unless ip\ — and / = 0. 



3 L p — L q estimates 

The representation formula (1.11) of Theorem 1.1 can be used to reproduce for the solutions of the wave 
equation in Einstein & de Sitter spacetime some important properties which possess the solutions of the 
wave equation in Minkowski spacetime. Among them there are estimates of the norm of solution in various 
functional spaces, such as L p , Sobolev spaces, Besov spaces and others. These estimates provide a useful 
tool to prove local and global in time existence theorems [24], [27], [34], [35]. 

In this short note we derive such estimates in the Lebesgue spaces only. First we remind these estimates. 
If n > 2, then for the solution v = v(x,t) of the Cauchy problem for the wave equation in Minkowski 
spacetime 

v tt -Av = 0, v{x,0) = <p(x), v t (x,0)=0, (3.1) 
with ip(x) e C ( 5 > °(]R™) one has (see, e.g., [3], [21]) the following so-called L p — L g decay estimate 

||(-A)-«(-»*)IU«(R») <Ct 2s - n ^-^y\\ LP{Rn) for all t > 0, (3.2) 

provided that s > 0, 1< p < 2, ± + \ = 1, and \{n + 1)(± - \) < 2s < n(± - \). 
Then, for the solution v — v(x,t) of the Cauchy problem for the wave equation 

Vtt - Av = 0, v(x,0)=0, v t {x, 0) = ip(x) , (3.3) 
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there is the LP — L q estimate 

||(-A)-' W (-,t)IU«(R»)<Ct 2a+1 - n( p-« ) ||v|Up(R») forall i>0, (3.4) 
under the conditions s > 0, 1< p < 2, ± + | = 1, and \{n + 1)(| - |) - 1 < 2s < n(| - |). 

The case of (ipo). According to Theorem 1.1, for the problem with ipi = and / = the function 
^ = ip(x,t) can be represented as follows: 

${x,t) = t^v^ix^t^-it-^tdtv^ix^t 1 ^). (3.5) 

Here for tp G Cq°(W"') the function w Vo (x, Si 1 / 3 ) coincides with the value v(x,3t 1 / 3 ) of the solution v(x,t) 
of the Cauchy problem (3.1). Hence for s > by means of application of (3.2) we obtain 

||(-A)-Xo(-,3t 1 / 3 )|U 9(R n ) < CtK 2s -«(|-|))||^ || iP(Rn) , 4 > o. 
To estimate the second term of (3.5) we apply (3.4) with s > : 

||(-A)- s (9 r ^ )(.,3t 1 / 3 )|| i , (R „ ) < Ct^ 2s+1 - n ^-^)\\ A Vo |Up(r»), * > 0, 

provided that ±(n + 1)(± - ±) - 1 < 2s < n(± - ±). Consequently if s > 0, 1 < p < 2, ± + ± = 1, and 
|(n + l)(i — i) < 2s < n(i — i), then for the problem with ^ = and / = we obtain 

||(-A)-'V(-,t)IU«(R») < Ci- 1+ K 2 -"(|-|))||^ || i!)(Mn) + Ct-i^( 25+1 -"^-|))|| A Vo ||ip(R») 

< C , *i( 2 '- n( p-i ) ) (t _1 || V o||iP(R»)+t~*|| A V0 ||iP(R»)), t>0. 

Thus, we have proven the following proposition. 

Proposition 3.1 Suppose that s > 0, 1 < p < 2, ± + ± = 1, and ±(n + 1)(± - ±) < 2s < n(± - ±). Tften 
i/ie solution tp — tp(x,t) to the problem 

j $ tt -t- 4 /3 A ^ + 2t-Vt = 0, O0,i6E", 

] lim ty(a;,i) = ¥>o(aO, lim (t-0*(M) + V^OM) +3i~ 1/3 A^(i)) = 0, xel", 
with tpo € Cq°(R") satisfies the following estimate 

||(-A)-V(-,t)IU.(R») < ^(^-"(l-D) (i-i||^o||L P ( K n) + ||A Vo |k I .(Mn)), t>0, (3.6) 
wit/i i/ie constant C independent of ip . 

The case of (<pi). For the problem with = and / = the function V = i>(x,t) due to Theorem 1.1 
can be represented as follows: 



ip(x,t) 



3 f 1 

- v Vl (x,4>(t)s){l - s 2 )ds, i£l",i>0. 
2 Jo 



Then we obtain for s, n, p, and g such that 2s — n(^ — ±) > — 1, the following estimate 
||(-A)-V(-,t)IU«(R-) < |jT IK-AJ-'^O.^H^^a-a 2 )^ 

< c AK^-d-i^s^-d-i^n^iu^a-s 2 )^ 

Jo 

< C^-^-i))|| Vl || iP(Rn) f s^-^-^il-s^ds. 

Jo 

Thus, in this case we have proven the following proposition. 
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Proposition 3.2 Suppose that s > 0, 1 < p < 2, i + ^ = 1, 2s - n(i - ^) > -1, and 5(71 + l)(i - ^) < 
2s < n(| — |). TTien f/ie solution ip — ip(x,t) to the problem 

J lira ttp(x,t) = 0, lira (tip t (x,t) + ip{x,t)) = ipi{x), x e R™, 
wii/i (£1 G Cg (K rl ) satisfies the following estimate 

||(-A)- s VM)||l,( R -) < CtK 2 -"(|-|))||^|| iP(Rn) , t> 0) (3.7) 
with the constant C independent of ipi . 

The case of (f). According to Theorem 1.1, for the problem with tpo = and ipi = the function 
^ = ip(x,t) can be represented as follows: 

o rl pl-b 1/3 

ip(x,t) = -t 2 dh dTbv f (x,3t^ 3 T;tb)(l + b^ 3 -t 2 ) . 

2 Jo Jo 

Consequently, for the problem with ip = 0, (pi = 0, and the function / satisfying conditions of the theorem, 
by application of (3.2) we obtain 

||(-A)->(-,i)|U< I (Mn ) 

< f 2 JdbJ dTb\\{-A)- s v } i-,it l /\-tb)\\ Lq(m {l + b 2 / 3 -T 2 ) 

,1 ,1-b 1 / 3 

I db / dr6tK 2 -«(M)) T 2 -«(i-|)||/(,t & )|| iP(R „ ) (l + & 2 /3_ T 2). 

Jo 



Ct 2 



For a = 2s — n{- — -)>—! one has 

,1-b 1 ' 3 



ri- — o 9 / \ a+1 / \ 

J ^(l + tV-Odr- » (!->■/») (l + , a + 



Hence, 



-A)-V(-,t)IU.(R») 

,1 ,1-b 1 / 3 



< Ct 2 +^-^-^) b\\f(;tb)\\ LP{Rn) db T 2S - n( ^\l + b 2 /*-T 2 )dT 

Jo Jo 

< C„, p , g , fl t 2+ K^-«(i-i )) jT 1 611/(^6)11^^) (l - 6V3) a+1 (1 + (a + i )6 V3 + 6 2/3) 

, g , s t 2 +K 2 -«(|-D) J' b\\f(;tb)\\ LP(Rn) (l - 



< C n , p . a J^*-'"7 j b\f(..lb)\\,, m , M -b> ) dh. 
Thus, in this case we have proven the following proposition. 

Proposition 3.3 Suppose that s > 0, 1 < p < 2, ± + ± = 1, 2s - n(± - |) > -1, and \{n + 1)(± - ±) < 
2s < n(- — -), and i/iai </ie function f satisfies conditions of Theorem 1.1. Then for the solution ip — tp(x, t) 



• 1 _ 1 
to ffte problem 



f iptt - t~ 4/3 A V + 2t~ Vt = /(a:, t), i > 0, 

] lim tr/>(ar, t) = 0, lim (tip t (x,t) + ip(x,t)) = 0, iel°, 

the following estimate 

||(-A)-'V(-,t)||L«(R») < C ms ti( 2s - n ^-^) / T\\f(;T)\\ LHRn) dT 

Jo 

holds with the constant C n , p , q ,s independent of f. 
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